APS/123-QED 



Finite Element Formulation of the Bloch Equations with Dipolar Field Effects 
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A Galerkin finite element (FEM) formulation for the Bloch equations with dipolar field is presented 
which makes possible the derivation of weak solutions to the Bloch equations. The FEM formulation 
has the advantage that the equations of motion are local in real space, eliminating the global 
truncation errors associated with calculations of the dipolar field in Fourier space. The dipolar field 
and other geometric parameters are calculated only once, before the simulation, and used as an 
initial condition rather than re-calculated at every time step of some numerical integration. 

PACS numbers: 76.60.Jx 



I. INTRODUCTION 

In many modern nuclear magnetic resonance (NMR) 
and imaging (MRI) experiments, the Bloch equations 
with dipolar field arc required to describe the observa- 
tion of subtle phenomena [J, 0, i, i, 0, 0, 0, 0, 0, EH El 
EH EH EH- F° r example, in EDM experiments, dynamic 
instabilities in highly polarized liquid xenon have been 
modeled in this manner [4|. Structural imaging at sub- 
voxel resolutions in MRI is possible using the long-range 
magnetic dipole field @, 0, 0, EH E3 ■ In all those cases, 
considerable insight in understanding the spin dynamics 
has been obtained from numerical simulations 0, 51, EH ■ 

The calculations generally involve the integration of 
the Bloch equations with a dipolar field. These are a 
set of nonlinear first-order partial differential equations 
that contain non-local integral operators in the calcula- 
tion of the derivative. The dipolar field integral intro- 
duces an 0(N 6 ) operation when calculated at each point 
in space, making the calculations very computationally 
intensive. For this purpose, Enss and Warren introduced 
a dramatic speed-up by calculation of the dipolar field 
in /c-space |15| . Any such simplification, tends to intro- 
duce global truncation errors into the numerical integra- 
tion, because the dipolar field must be calculated at every 
time step using solutions that contain errors accumulated 
from each of the previous steps. The cumulation of global 
truncation may explain the inability to obtain numerical 
solutions at long evolution times. 

We introduce a finite clement (FEM) formulation 
which may alleviate these problems, because the dipo- 
lar field need not be recalculated at each time step. This 
Galerkin formulation of the Bloch equations with dipolar 
field was first presented in reference [l6| . 



II. BLOCH EQUATIONS WITH DIPOLAR 
FIELD 

The secular part of the dipolar field is [H EH : 
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Although we use the secular approximation, which is 
valid in high magnetic fields, the calculation herein may 
also be done for the low field case using the complete 
dipolar field including all non-secular terms [111, EH : 
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The exact form of the dipolar field is unimportant here. 
Chemical shift offsets or radiation damping effects may 
also be readily added. 

In addition to DDF, diffusion and relaxation, flow ef- 
fects can also be accounted for (in this section we write 
M = M(r, t) and B = B d (r, t) for the dipolar field): 
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This is the set of partial differential equations that is 
normally integrated by quadratures to provide numerical 
solutions 0, EH- 

The diffusion term V 2 M requires the existence of the 
second derivative. This is a rather strong requirement 
which is not usually needed. In the theory of heat con- 
duction, weak solutions to the heat diffusion equation 
dT/dt = kV 2 T often model physical situations well. In 
the next section, we investigate the existence of weak so- 
lutions. 
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A. Weak solutions 



Let C M 3 be the diffusion region and dfl be its 
boundary. The Bloch equations with dipolar field B(r) 
and magnetic field offset 5B(r) (to allow for magnetic 
field inhomogeneities) are: 



dM - -5 5 + M 2 e 2 

M - M 3 9 -j. , s 

+ — °— ±e 3 + D\7 2 M (3) 

J-i 

are defined for r e and supplemented by appropriate 
boundary conditions on dQ, for example n«VM = on 
dfi in the case of reflective boundaries. As usual, B is 
the dipolar field. In component form (component i) this 
reads: 
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Let V be the vector space of test functions on fi, i.e. 
all functions that are continuous, satisfy the boundary 
conditions on dfl and whose first derivative is piecewise 
continuous on Q (see Quarteroni [l7| for details). Take 
a test function Vi(r) £ V, multiply Eq. (TJJ by Vi and 
integrate over f2: 
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We proceed to eliminate the second derivative in the 
diffusion term. This is done by using the vector identity 
V-(/sO = (V/H + /(V.ff) with / = Vi and g = VM ; . 
The diffusion term becomes: 
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where in the second line we made use of the divergence 
theorem. The surface integral vanishes by virtue of the 



boundary conditions n«VAf = on dil. Thus, the diffu- 
sion term reduces to 



D I v t V 2 M t d 3 r = ~D I (VMi-Vv t )d 3 r (7) 
Jn Jn 

We have derived the following useful result: 

Definition II. 1 (Weak form of the Bloch equations). 
The weak form of the Bloch equations — with distant dipo- 
lar field, arbitrary distributions of resonance frequency 
offsets, and diffusion within reflective boundaries — is 
given by (no summation on i): 

J ^jf v i d3 r = 7 ]T e ijk J Mj(B k + 6B k ) Vl d 3 r 
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where Vi is any test function. Any function Mi which 
satisfies this equation is called a weak solution. 

We note that the degree of differentiability of Mi has 
been reduced by half. For the purpose of solving this 
problem numerically, the Galerkin method can be intro- 
duced. 

Definition II. 2 (Galerkin formulation). Let Vh be a fi- 
nite dimensional subspace of V of dimension Nh . The 
Galerkin finite element method consists of finding three 
functions u^(r, t) £ Vh (i = 1, 2, 3) such that: 
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(There is no summation on the index i.) The Galerkin 
formulation is the matrix form of this problem. 

Remarks. 

1. The Galerkin method is only an approximation be- 
cause the solution is expanded in terms of a finite- 
dimensional basis of functions. 

2. The are the best approximations to Mi (in the 
L 2 sense) that the finite-dimensional vector space 
Vh allows, i.e. u\ is the orthogonal projection of 
Mi on the subspace Vh- 
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Next, we pick a set of basis functions {(fi} for the vec- In this basis, the dipolar field B k approximates to: 
tor space Vh and expand the solution in terms of this 
basis: [H| 



(10) 



where the functions Wi„(t) are the unknown coefficients 
we must solve for. Taking = ip m yields: 
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where a k = — 1 if k = 1, 2 or a k = 2 when fc = 3. Let- 
ting ||A nm || be the matrix whose elements are A„ m = 
J„ tPntPmd ' r, and assuming it is non-singular, we denote 
the inverse of ||A„ m || by ||C nTO ||. Substitution of -Bfc and 
multiplication of the expression by the matrix ||C nrn || al- 
lows us to solve for the time derivatives: 
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Then, with the following abbreviations, 
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the expression takes a particularly simple form, 



This is a set of coupled ODEs, involving linear and 
bilinear forms in the unknowns u^-, with constant coef- 
ficients R k n7m Cprm A nvai Htm G nra and T knra q. These 
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constants can be calculated a priori using the choice of 
basis functions {tfi}. One may then solve for the time- 
dependent coefficients Wi P (there are 3 x Nh of them) us- 
ing conventional ODE methods (see, for example, Quar- 
teroni fPf\). 

III. CONCLUSION 

We have shown how to derive the Galcrkin formula- 
tion of the Bloch equation. The main advantage to this 
approach is in calculating time-evolution of magnetic mo- 



ments in the presence of a strong dipolar field, in which 
case, we believe that error propagation from repeated 
dipolar field calculations should be mitigated [!, 0, [l(| • 
In contrast with the direct integration of the Bloch equa- 
tions using the Fourier method [3, [l5[ , there is no need 
to calculate a dipolar field at every point in time; hence, 
there is no propagation of global truncation errors associ- 
ated with the Fourier approximation of the dipolar field. 
This set of equations, while still non-linear, docs not in- 
volve a non-local integral operator in the calculation of 
the derivative at each time step. 
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